Replicate Periodic Windows in the Parameter Space of Driven Oscillators 
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In the bi-dimensional parameter space of driven oscillators, shrimp-shaped periodic windows are 
immersed in chaotic regions. For two of these oscillators, namely, Duffing and Josephson junction, we 
show that a weak harmonic perturbation replicates these periodic windows giving rise to parameter 
regions correspondent to periodic orbits. The new windows are composed of parameters whose 
periodic orbits have periodicity and pattern similar to stable and unstable periodic orbits already 
existent for the unperturbed oscillator. These features indicate that the reported replicate periodic 
windows are associated with chaos control of the considered oscillators. 
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I. INTRODUCTION 

The parameters of deterministic dynamical systems 
play an important role to specify the transitions between 
chaotic and periodic behavior. This parameter influence 
on the attractor transition can be represented in the pa- 
rameter space In particular, for both discrete- time 
[1, 6] and continuous dynamical systems 0, it is known 
that the set of parameters for which a system exhibit pe- 
riodic behavior are periodic windows immersed in chaotic 
regions, in the form of shrimps, in the bi-dimensional pa- 
rameter space. In the past two decades, periodic windows 
have been numerically obtained for a large number of ap- 
plied dynamical systems like lasers |8|. electronic cir- 
cuits |9l-[llj , mechanical oscillators [12[ , and also in pop- 
ulation dynamics [l3j |. Periodic windows have also been 
identified in experiments with electronic circuits [l4]-[l6j . 

Additionally, the bi-dimensional parameter diagram is 
an important tool to analyze the outcome of techniques 
used to control chaos for allowing a global overview of pe- 
riodic and chaotic behavior [17j. In literature, different 
approaches have been proposed to control chaos. Essen- 
tially, two main methods have been successfully employed 
in several dynamical systems, namely the feedback and 
the nonfeedback methods. The feedback methods are 
applied to maintain the trajectory in a desired unstable 
periodic orbit embedded in the chaotic attractor [l8|. In 
contrast, the nonfeedback methods eliminate the chaotic 
behavior by a slightly modification in the system dynam- 
ics [l9| . The suppression of chaos by applying an external 
weak perturbation is an interesting application of a non- 
feedback method. 

Suppression of chaos by adding a weak harmonic per- 
turbation has been numerically and theoretically re- 
ported specially for systems with an original harmonical 
driven [20j. For example, a weak harmonic perturbation 
has been used to suppress chaos in the forced Duffing 
oscillator [21], and a similar perturbation has been ap- 
plied to control chaos in a Josephson junction oscillator 
|22| . Considering that in |2CH22] the controlling of these 
two well-known systems were essentially accomplished for 



limited ranges of system parameters, a further parameter 
space analysis considering a weak harmonic perturbation 
is required. 

The main motivation of this article is to determine the 
alterations in the parameter space due to weak harmonic 
perturbations applied to continuous-time dynamical sys- 
tems. To investigate that we analyze how such pertur- 
bations modify the attractors and the parameter space 
of the harmonically driven Duffing and Josephson oscil- 
lators. For the first time in the literature, we observe 
that this harmonic perturbation replicates the periodic 
windows in the parameter space of these systems. More- 
over, we find evidences that new periodic orbits, whose 
parameters are in the new periodic windows, are similar 
to unstable periodic orbits embedded in the unperturbed 
chaotic attractor and periodic orbits already existing. We 
discover that the reported replicate periodic windows are 
associated with chaos control of the considered oscilla- 
tors. 

This letter is organized as follows: In Section [III we 
obtain the parameter space of the unperturbed Duffing 
oscillator and we establish the suitable parameters to im- 
plement the weak perturbation. In Section IIII( for the 
Duffing oscillator, we investigate the periodicity, shape 
and Lyapunov exponents of the perturbed periodic and of 
the unperturbed chaotic orbit. In Section HVl we present, 
for the Josephson oscillator, another example of replicate 
periodic windows caused by a weak periodic perturba- 
tion. Finally, in Section 5 we summarize our main con- 
clusions. 



II. PARAMETER SPACE STRUCTURE OF THE 
DUFFING OSCILLATOR 

The Duffing oscillator is a well-known model to de- 
scribe oscillations of a mass obeying a fourth order sym- 
metric potential [23|. This system has a large number of 
applications, specially to model physical systems, and has 
been extensively studied in theoretical, numerical and 
experimental approaches (23-[26|. Here, we consider a 
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simple version of the Duffing system, which describe os- 
cillations of a single-well potential. 

The time evolution of this system is determined by 
solution of the following dimensionless equation: 

x + cx + x 3 = /3 cos(cjt). (1) 

Here, the parameter c is the amplitude of the system 
damping, f3 is the forcing amplitude, and uj is the natural 
system frequency settled at lj = 1.0. We numerically ob- 
tain solutions of the Duffing equation by using a fourth- 
order Runge-Kutta method with fixed step h = 0.001. 

Additionally, to investigate transitions between chaotic 
and periodic attractors of the Duffing oscillator, the 
largest non zero Lyapunov exponent is obtained [27j for 
each point of a bi-dimensional grid of the system param- 
eters (parameter space). 

In Figure [TJ we show the bi-dimensional parameter 
space (c x 3) obtained for the Duffing oscillator (also 
shown in |25j). The Lyapunov exponent features of each 
point in the grid are represented by assigning different 
colors. In Figure [TJ blue color represents parameter sets 
for which the attractors are chaotic (positive Lyapunov 
exponent), while grey-scale represents parameter sets for 
which the attractors are periodic (negative Lyapunov ex- 
ponent). In black we represent points of the superstable 
lines corresponding to periodic attractors for which the 
largest non zero Lyapunov exponents is a minimun inside 
a periodic window. 

We observe, in Figure [TJ that the parameter sets for 
which the Duffing oscillator behaves periodically are in 
aggregated periodic windows (gray scale area). Those 
periodical structures, also known as shrimps, have been 
well described in literature p]-[3|. In Figure [2(b), the 
red circle and the black plus symbol indicate examples 
of parameter sets whose orbits are, respectively, periodic 
and chaotic; these two orbits will be further considered 
in Section [Till 

To show the different behavior of the two attractors 
associated with the two parameter sets marked in Fig- 
ure DIb), we obtain a stroboscopic map by collecting the 
velocity and the displacement at (Time-27r/o;). In Fig- 
ure Efa), we show the periodic orbit correspondent to the 
parameters represented in red inside the periodic window 
amplified in Figure [2(b). Similarly, in Figure [2fb), we 
show the chaotic attractor for parameters represented by 
a plus symbol in Figure [2(b). 

We call the attention, in Figure [2j that both periodic 
and chaotic attractors are in the same phase space region. 
Next, in Section HIH we analyze how these two orbits are 
modified by a weak periodic perturbation. 

III. INTRODUCING A WEAK PERTURBATION 
IN THE DUFFING OSCILLATOR 

We introduce a weak harmonic perturbation by adding 
a second harmonic term in the original Duffing driven. 



The perturbation amplitude is taken as a control param- 
eter of the system. On the other hand, the perturbation 
frequency is settled in an integer ratio of the original 
system frequency. Similar procedures have been already 
reported in literature [20|-[22j. So, the time evolution of 
this perturbed system is determined by solution of the 
following dimensionless equation: 

x + cx + x 3 = /3 cos(cjt) + asin(Qt). (2) 

Here, a is the weak perturbation amplitude (a << j3) 
and ft is the perturbation frequency fixed at ft = 2oo. 
Other rational multiples of uo could be used to Vt [I?], Hl| . 

In Figure O we present the parameter space of the 
perturbed Duffing oscillator for two different perturba- 
tion amplitudes a = 0.04 and a = 0.08. By comparing 
the parameter spaces of Figure [2(b) and Figures [3](a,b), 
we note the existence of replicated periodic windows. In 
particular, the central periodic window of Figure QJb) 
is duplicated in Figure [3l Moreover, the new periodic 
windows appear slightly displaced in Figures 02(a) and 
02(b). Thus, in Figure 02(b), for the perturbation ampli- 
tude a = 0.08, the black plus marked point is inside the 
replicated periodic window. In other words, the chaotic 
attractor showed in Figure [2fa) has been changed to a 
periodic one. 

Next, we discuss the alteration produced by the weak 
perturbation on the chaotic orbit shown in Figure [2fb). 
For this purpose, in FigureH^a), we show the stroboscopic 
phase space of the new periodic orbit (in black plus sym- 
bols), that substitutes the former chaotic orbit shown in 
Figure Efb). Additionally, to compare similarities be- 
tween perturbed and unperturbed orbits, in Figure Ufa), 
we present (in red points) the unperturbed periodic or- 
bit whose parameters are indicated in red in Figure QJb). 
In Figure 2Kb), we show the convergence of the largest 
non zero Lyapunov exponent of the two orbits of Fig- 
ure Ufa). We recognize, in Figure E](a), that the new 
perturbed periodic orbit and the previous periodic orbit 
have the same periodicity and similar pattern. Moreover, 
the largest non zero Lyapunov exponent of both orbits 
also have a similar convergence. 

We investigate the relation between the new perturbed 
periodic orbit and a possible unstable periodic orbit em- 
bedded in the chaotic sea. Thus, for the parameters 
marked by a plus symbol in Figure [2(b), we integrate 
the unperturbed equation for a short integration time, 
starting with initial conditions settled in the unperturbed 
periodic orbit with parameters marked by red circle in 
Figure [2(b). 

In Figure the blue triangle symbol denotes the 
chaotic orbit in the short integration time, the red circle 
symbol denotes unperturbed existent periodic orbit, and 
the black plus symbol denotes the perturbed periodic or- 
bit. In Figure we observe that for a short integration 
time, the chaotic orbit remains in the neighbor of the 
perturbed and the unperturbed periodic orbit. 
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FIG. 1: (Color online) (a) Bi-dimensional (c x f3) parameter space of the unperturbed Duffing oscillator, for uj = 1.0. (b) 
Magnification of the squared area. The red circle indicates periodic behavior at f3 = 10.6780 and c = 0.2670, the black plus 
symbol indicates chaotic behavior at /3 = 10.7020 and c — 0.2670. 
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FIG. 2: (Color online) Stroboscopic phase space of the Duff- 
ing oscillator, (a) Periodic orbit correspondent to the pa- 
rameters f3 = 10.6780 and c = 0.2670. (red circle in Figure 
1(b)). (b) Chaotic attractor correspondent to the parameters 
f3 = 10.7020 and c = 0.2670 (black plus symbol in Figure 
1(b)). 



The analysis of results shown in Figures |4] and [5] in- 
dicate that the periodic window replication described in 
this work gives rise to new periodic windows whose pa- 
rameters correspond to both stable and unstable peri- 
odic orbits with the same periodicity and pattern found 
in previous existing unperturbed oscillator. 



IV. WEAK PERIODIC PERTURBATION IN 
THE JOSEPHSON JUNCTION OSCILLATOR 

The superconducting Josephson junction is usually 
modeled by an electronic circuit equation [28]. Consid- 
erable efforts has been devoted to understand the onset 
of chaos in thi s sy stem and its critical parameters for 
this transition |29i-l31|. The system is described by a 
pendulum-like dimensionless equation (22j : 

4> + G(j) + sin (j) = I + A sin(cjt). (3) 

Here, the parameter G gives the amplitude of the 
system damping, / is the direct current component of 
the circuit, uj is the natural system frequency settled at 
uo = 0.25, while A is the alternating current component 
of the circuit and the system forcing amplitude. 

In Figure [6j we obtain the bi-dimensional parameter 
space (G x A) of the Josephson oscillator given by Equa- 
tion [3j In this figure, blue points correspond to positive 
Lyapunov exponent (chaotic behavior), while grey- scale 
points correspond to negative Lyapunov exponent (pe- 
riodic behavior). The black region corresponds to the 
superstable lines. 

We verify, in Figure [6j the existence of aggregated pe- 
riodic windows in the parameter space of the Josephson 
oscillator. As in Figure DJb), the red circle and the black 
plus symbol indicate in Figure [6fb) examples of parame- 
ters for which the correspondent orbits are, respectively, 
periodic and chaotic. 





10.90 



FIG. 3: (Color online) Bi-dimensional (c x /3) parameter space of the perturbed Duffing oscillator for uj = 1.0 and Q = 2.0. 
The perturbation amplitude is a = 0.04 in (a) and a — 0.08 in (b). 
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FIG. 4: (Color online) (a) The red circle symbol denotes the 
periodic orbit for /3 = 10.6780, c = 0.2670 and a = 0.0, the 
same orbit of Figure 2(b). The black plus symbol denotes 
the controlled orbit for f3 = 10.7020, c = 0.2670, Q = 2.0, 
and a = 0.08. (b) The red line indicates the largest Lya- 
punov exponent of the unperturbed periodic orbit. The black 
line indicates the largest Lyapunov exponent of the controlled 
orbit. 
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FIG. 5: (Color online) The blue triangle indicates the chaotic 
orbit under a short integration time for /3 — 10.7020, c — 
0.2670 and a = 0.0. The black plus symbol indicates the 
controlled periodic orbit for f3 = 10.7020, c = 0.2670, a = 0.08 
and Q = 2.0. The red circle symbol indicates the unperturbed 
periodic orbit for = 10.6780, c = 0.2670 and a = 0.0. 



term is added to the original one: 



In Figures [T^a) and [71(b) we present, respectively, the 
periodic and the chaotic orbits for the two parameter sets 
indicated in Figure [6fb), to be further analyzed in this 
section. 

For the Josephson oscillator, similarly to the perturbed 
Duffing oscillator (Considered in Section 3), a harmonic 



<f> + G<j) + sin cj) = I + A sm(ut) + a sin(fit) . (4) 

where, a is the perturbation amplitude and Q is the 
perturbation frequency fixed at Q = uj/2. 

In Figure [8] we show the parameter space of the per- 
turbed Josephson equation. Similarly to what was no- 
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FIG. 6: (Color online) (a) Bi-dimensional (G x A) parameter space of the unperturbed Josephson oscillator, the natural 
system frequency is settled for uj = 0.25, the direct current component is settled for / = 0.905. (b) Magnification of the squared 
area. The red circle indicates the parameters A = 0.5289 and G — 0.4403, the black plus symbol indicates the parameters 
A = 0.5488 and G = 0.4416. 
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FIG. 7: Stroboscopic phase space of Josephson oscillator, (a) 
Periodic orbit correspondent to the parameters A = 0.5488 
and G = 0.4416. (red circle in Figure 6(b)). (b) Chaotic 
attractor correspondent to the parameters A = 0.5488 and 
G = 0.4416 (black plus symbol in Figure 6(b)). 



turbed oscillator. 

Next, to show how these orbits are modified by the con- 
sidered perturbation with amplitude a = 0.01, we show 
in Figure [9] the stroboscopic phase space of the perturbed 
Josephson oscillator with three orbits. Thus, in Figure [9j 
the black plus symbol represents the controlled periodic 
orbit indicated by a black plus symbol in the parameter 
space shown in Figure [5(b), with the same parameters 
of the former chaotic orbit marked with the black plus 
symbol in Figure [6(a). The blue triangles indicate the 
chaotic orbit (whose parameter set is indicated in Fig- 
ure Eft))) under a short integration time for the initial 
condition fixed on the periodic attractor. The red circle 
symbol denotes unperturbed periodic orbit. 

Thus, the results shown in Figures [6] - [9] confirm the 
evidence, found for the Duffing oscillator, that the new 
periodic orbits are similar to previous stable and unstable 
periodic orbits of the unperturbed oscillator. 



V. CONCLUSIONS 



ticed in Figure [3] for the Duffing oscillator, in this case 
the existence of replicate periodic windows in the mod- 
ified parameter space is also observed. Moreover, com- 
paring Figures [6] and [5J we see that the parameter set 
indicated by a plus symbol in these two figures represent 
a chaotic attractor in Figure [6(b), for the unperturbed 
oscillator, and a periodic one, in Figure [5(b), for the per- 



We investigate the control of chaos for two driven os- 
cillators by a weak harmonic forcing. To identify peri- 
odic and chaotic regions in the bi-dimensional parameter 
space, we compute the largest non zero Lyapunov ex- 
ponents for the attractors in the considered parameter 
ranges. We identify shrimp-shaped periodic windows im- 
mersed into a chaotic region. The parameter space is 
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FIG. 8: (Color online) Bi-dimensional (G x A) parameter space of the perturbed Josephson oscillator, the natural system 
frequency and the perturbation frequency are settled, respectively, for uj = 0.25 and Q = 0.125. (a) The perturbation amplitude 
is settled for a = 0.005. (b) The perturbation amplitude is settled for a — 0.01. 




FIG. 9: (Color online) The blue triangle indicates the chaotic 
orbit under a short integration time for A = 0.5488, G — 
0.4416 and uj = 0.25. The black plus symbol indicates the 
controlled periodic orbit for A = 0.5488, G = 0.4416, uj = 0.25 
and Q = 0.125. The red circle symbol indicates the previous 
periodic orbit for A = 0.5488, G = 0.4416 and uj = 0.25. 



much modified whenever a weak amplitude forcing is ap- 
plied. New similar periodic windows arise in the neigh- 
borhood of the original windows. We verify that peri- 
odic orbits are similar (and with the same periodicity) 
for parameters inside the original and the new periodic 
windows. These results are similar to those obtained for 
the driven impact oscillator [l7| . 

By analyzing stroboscopic maps of unperturbed and 
perturbed attractors we find evidences that the new re- 
ported periodic windows are formed by parameters for 
which the observed new periodic attractors are similar 
to preexisting stable and unstable periodic orbits. There- 
fore, we conjecture that the replicate periodic windows 
reported in this work are associated with chaos control 
and reproduce further other periodic orbits of the con- 
sidered oscillators. 
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